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1.  Introduction 

Successful  prediction  of  the  amount  of  error  incurred  in 
gunner  tracking  of  a  moving  target  has  become  an  important  phase  in  the 
development  of  weapon  systems.  A  great  deal  of  data  in  engineering, 
economics,  and  the  natural  sciences  occurs  in  the  form  of  time  series 
where  observaciy.ns  are  dependent  and:  where'  the  nature  of  this  dependence 
is  of  utmost -importance . 

The  purpose  of  this  report  is  to  present  an  evaluation  of  the 
gunner  error  data  described  in  Section  2?  thereby  defining,  a  time  series 
model.  This  modei  and  its  parameter.^  can  be  used  to  simulate  data  for 
future  problems  of  a  similar  nature  of  may  be  used  as  a  subroutine  to 
-missile  flight  simulation. 

2.  Data  Description 

Tests  to  determine  the  gunner  tracking  error  characteristics 
were  conducted  at  Redstone  Arsenal  during  the  period  13  through  18, 

July  1972.  The  King,  Air ,  a  twin  engine  Beechcraft,  was  the  target, 
utilized  for  these  tests. 

A  16mm  film  camera  was  attached  to  the  monocular  output  of  the 
tracker  unit.  This  output  presents,  the  same  view  to  the  film  camera  as 
,che  binocular  output  .presents  to  the  gunner. 

There  were  four  gunners  who  participated*  in  the  tests.  They  were 
.instructed  to  track 'the  centroid  to  the  target  aircraft  when;  details 
were  riot  resolvable..  Wheri*  resolvable  they  were- to  track  the  inter¬ 
section  at  the  wing,  and  fuselage.  The  amount  of  error  was  shown  to  be 
independent  of  individual  gunner,  that  is,  there  was  no  statistical 
significance . 


3.  Model  Building 

This  section  discusses  the  time  series  model  building  for  the 
gunners:'  error  data.  After  examining  all  the  data  available,,  we  conclude 
that  the  data  forms  a  stationary  time  series  except  at  the  beginning 
where  a  transient  occurs  during  acquisition,  and’  at  the  end  where  a 
transient  is  introduced  by  the  simulated  missile  in  flight  signal. 

Runs  with  too  few  data- were  eliminated.  The  total  number  of  runs  was 
then  143.  A  few  nonstationary  data  can  also  be  seen.  They  occupy 
13.29  percent  of  the  'total. 

When  the  data,  are  recorded  with  equally  spaced  time  intervals,  we 
generally  use  a  linear  time  series  model  to  fit  the  data.  A  commonly 
used  model  for  univariate  time  series  can  be  written  as 
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(1.1) 


\  “  **  =  *i^t-r  ”  +  ^2at-2  "  ^  +  *••• +  %(Y't-P  "  fi) 


^  ^ .  •  0 1 3 .  t-|  —  *  *4.  —  Q-  3 

t  1  t-1  q  t-q 


where 


subscript  t  =  time 

Y  ■?  the  value  of  the  time,  series  at  time  t 

p  =  the  expected  value  of  Yfc 

afc  =  a  white  noise  process,  i.e.„  a:  is  independent, 

2 

identically  distributed  N(0,  a  ) 

3 

p,  q  =  two  parameters  depending  on  the  properties  of  a 
particular  -time  series. 

Model  (1.1)  is  called  a  mixed  model  with  autoregressive  and  moving 
average  components.  It  has  beenwidely  used  in  practice  with  fruitful 
results  (see  e.g.,  Box  and  Jenkins  fl),  Fuller  and  Tsokos  [2];,  Cleveland 
[3,  4]  ,  and  Box  et.  al  [5]),.  The  intuitive  idea  behind  the  model  (1.1) 
is  the  assumption  that  the  present  value  Y  depends  on.  the  values  of 

Yt  the  .near  past,  i.e< j  Y^,  'Yfc_2»  Yt-p’  ^his  is  the,  autore¬ 

gressive  component 

(Yt  -  tx)  =  -  J)  +  <M¥t_2  "  *»)'  +  ...  +  Op:(Y^p  -  u)  . 

The  moving  average  component  a.  -  6, afc  r  . ..  -  0  a-  indicates  that 

c  l  t-i  q  t-q 

’the  present  value  Y  depends  not  only,  on  the  present  noise  a  ,  but  also 

c  t 

the  previous  noise  at_^,.  ...,  a-^  .  This  is  reasonable  since  the  noise 

will  not  diminish  very  rapidly  in  real  situations.  The  noise  prolongs 
its  influence  on  Yt  for  a  certain  period. 

In  practice  when  time  series  data  are  given,  a  model  of  the  form 
(l.l)  can  generally  'lie  built.  The  detailed  procedure  has  been  given  in 
Box  a,nd  Jenkins  [1],  There  are  four  main  steps. 

a •  Model  Identification 

In  this  first  step,  autocorrelation  coefficients,  partial 
autocorrelation  coefficients,  and  inverse  correlation  coefficients 
(e.g.,  Cleveland  { 3 J )  are  used  to  determine  the  values  of  p  and  q  in 
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model  (1.1) .  The  value  p  is  called;  the  order  of  the  autoregressive 
component  and  the  value  q  is  called  the  order  of  the  moving  average, 
component  in  a  mixed  model  (1.1): 

b .  Parameter  Estimation, 

After  the  values  of  p  and  q  have  .been  determined,  there 

are  p  +  q  +  2  parameters:  u,  4>.,  <!>„  ...,  ,  6,,  .  ..,,-,0  and  , the 

2  ^  *  ^p  l  ■  q 

variance  a  of  to  be  determined.  The-  method  used  to  estimate  e’s 

'  a  t  ,  "  " 

and  0's  has  been  described  in  Box  and  Jenkins'  (Chapter  7,  [11), 
Clevensoh  [6] ,  and  Parzen  [ 7] .  The  main  technique  is  the,  maximum 
likelihood'  estimation.  Generally,  the  calculation  needs  the  Help  of 
spectral  density  estimation  [7]  or  nonlinear  least  squares- estimation 


2 

The  estimated',  values  Q.  o, and  a  of  the  parameters 

TV  **  /Vfl  '  *  '  ’ 

9  * 

14,  o»  6 a  and.  o~  j.  respectively  ,  are  not  generally  equal  .to.  .the- >rea  1  -value 
of  these  parameters.  The  model  with  estimated  parameters 


(Yt  ~ =  %(Yt-i  “ oy  +  •  • +  V(Yt-p  * a)  +  a.t  ”  0Vt-i 


....  -  §  a, 

q  t-q 


11.2) 


may  not  fit  dihe  original  data  well..  Diagnostic  checking  determines 
whether  our  estimated  model  fits  the  data  well.  The  residual  process- 
[SjJ.  is  examined.  If  the>  [£  ,).  is  close  to  a  white  noise  .process,  the 

model  is  considered  to  be  adequate  and  the  whole  model  building  procedure 
is  over .  Otherwise,  we'  go  -to  the  next  step; 

d .  Modification  of  the  Model 

If  the  model  we  built  is  found'  inadequate  through  the 
diagnostic  checking,  we  will  try  to"fit  the  data  by  a  new  modified 
model.  Generally,  the  residual  process  •[£  ]3  will  reveal  some  information 

on  how  the -model  should  be  rebuilt.  In  most  cases,,  a>  pair  ,of  new  values 
of  p  and  q  will  be  obtained.  Using  these  new  values  of  p  and  q,  we 
undergo  steps  b.,  c.,  arid  d.  for  this  new- model  building. 

All  the  four  steps  have  been  carefully  followed  for  building  the 
gunners'  error  data  model.  For  the  (apparently)  stationary  time  series, 
with  azimuth  and  elevation  both  cor.nted,  the  total  number  of  realizations 
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/  ■ 

was  248.  Each  time  series  of  azimuth  and'  elevation  is  run  separately 
(•Tables  1  and  2) .  Sixty- two  percent  of  the  stationary  series  can  be 
fitted  well  by  a  third  order  autoregressive  process  [p  =  3,  q  =  0  in 
model  (1.1)],  i.e,.,  .  - 

Yt  '  »  =  ♦r(Yt-l  "  ^  4  *;2(Yt-2  ‘  +  4,3(Yt-3  "  ^  +  at  ' 

(1.3) 

A  few  data  -can  not  be  fitted  welt  by  (4.3);.  they  are  fitted  by  a  more 
complicated  model r  These  models  and  their  percentages  of  the  total 
data,  are  given  in  Table  4.  Due  to  the  biological  and  psychological 
differences  among  gunners,  there  are  variations  in  these  parameters. 

The  means  and  variances  of  these  parameters^  are  also  given  in  Table  1. 

TABLE  1.  GUNNER'S  EFROR  MODEL  FOR  AZIMUTH 


/ 

A 


V 


TABLE  2.  GUNNER'  S  ERROR  MODEL  FOR  ELEVATION’  , 


General  model  (3rd  order  autoregressive  process)  (60.48%). 


a 

-  v 

*2  *3 

2 

°a 

Mean  -0.0515 

0.3692 

0.2165  0.1*448 

0.0045 

Variance  0.0124 

0.0188 

0.0051  0.0057 

0,0001 

'Special  Model 

Mean 

Variance 

1)  $4  *■  o 

(15.32%), 

0.1535; 

0.0014 

2)  <t> 5  *  0 

(8.87%) 

0.1255 

0.0088 

3)  0 

(7.25%) 

0.1125 

0.0039 

4)  o  7  *  .0 

(4.84%) 

0.1221 

0.0087 

5)  c8  /.  0 

(0 .81%) 

0.1497 

0 

6)  09  *•  0 

(0.81%) 

0.1294 

0 

7)  .,10  ,  0 

(0.81%) 

0.1178 

P 

8  )  0  ^  r  0 

(0.81%)' 

0.0573 

0 

A  question  arises  whether  the  azimuth >error  and-  elevation  error 
are  dependent  on  each  other  during  a.  gunner's  aiming.  The  data  show 
that  we  can  consider  the  azimuth  error  and  elevation  error  to  be  twc 
independent  -processes .  The  following,  procedure  is  followed i 


A  general  model  describing  the  relation  between  two  time  series  is 

a  linear  transfer  function  model..  Let  X_  be  the  time  series  of  azimuth. 

t 

and  Y_  be  the  time  series  of  elevation.  A  linear  transfer  function 

t  ..... 

model  can  be  written  as 


,(Yt  “  V  =  «l<Yt-l  -  V  +  •••  +  <VYt-m  -  V 

f  MXt-l  “  ***>  *  +  ^n(Xt-n  ~  |Ix)  +  Nt 


(1.4) 
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•IU.P'* U.UPW 


WWWPBi 


where 


V=E(Yt) 

^  =  E(Xt) 

N-  =  a  noise  process 


m,,  n  =  the  numbers  of  past  values  of  X  and  Y  on  which  the  present 

XT.  J:. _ I  „  \  %  t  ' 


Y;  depends. 


Intutitively ,  model'  (i. 4)  indicates  that  the  present  azimuth  value  Y 

may  depend,  on  the  previous)  values  of  both  azimuth  and  elevation.  This 
model  has  been  used  in  many  practical'  situations  and  gives  good'  resuits 
(see.  e.g.,  Box  and  Jenkins  Since  we- ’have  already  found*  a  good 

model  for  Y^  in  the  previou's  model-  buildings.,  we  may  combine  the,  Y 

model  and  (1.4‘)  and  have 


0-5) 


where  a-  is  the  noise',  .process  from  the  model  of  -Y<  ,  Since  a  is  a 

white  noise,  process,  the  values  j3's  can  be  easily  estimated  (Box  and 
Jenkins  (ij  p .  380) . 

•'n-  attempt  has  been 'made*  to  fit  all  the  .corresponding  pairs -of 
azimuth,  error  data  and  elevation  error  data  by  model  (1.4).  Except  for 
a  few  exceptions  (1  perce  vt  of-  the  total)*,  the  6  values  are  very  small 
(dess  than.  0 .05  for  all  (3^,.  §23)  •  Hence,  consider  that  the 

error  in  elevation  has  no  significant  influence  on  ithat  in  azimuth.  A 
similar  model  fitting  by  replacing  X  by  Y  and  Y  by  X  in  (1.4)  has  also 
been  run  for  ail  pairs  <of  data*.  An  independence  relation  is  also 
obtained  here-.  Hence,,  we  conclude  that  there  is  <ho  significant 
dependence  between  azimuth  error  and  elevation  error.. 


4.  Simulation  Procedure 

In  order  to  simulate  the  total  performance  of  a  guided  missile 
system  with  a  man  in  the.  loop,  we  may  use  the  gunner's  model  described 
in  the  previous,  section.  Considering  the  nonrepeatability,  of  man's 
reactions,  it  must  be  realized  that  for  any  single  simulation  the  error 
model  will  not  give  the  same  results  as  given  by  man.  However,  man's 
behavior  oh  the  average  should  a,gree  with  that  of  the  error  mode. 
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Simulation -of  a  gunner’s  behav.ior  may  be  performed  as  follows: 

a)  Choose  2  random  numbers'^  and  y 2  in:[0,  1]  .  is  used  to 
construct  azimuth,  error,  if 

7^e  [0,  0.6290].,  a  -third  order  autoregressive  model  will  be 
used, 

7^,c  [0.629.1,,  0.8064],  a  fourth  order  autoregressive  model  with 

,0^  t-  0  will  be  used5, 

7^e  [0.8065,  0.8548] j  a,  fifth  order  autoregressive  model  with 

0,-  *  0  will  be  used, 

7^fc‘-J0.8549,  0.9113],,.  a  sixth  order  autoregressive  model  with 

4>g  vi  0  will,  be  used-, 

7,,e  [0.9li4,  0.9435].,  a'  seventh,  order  autoregressive  model 

with  -/  0  will  be  used, 

7^c  [ 0. 9436,,  0.9516] ,  an  eighth  order  autoregressive  model  with 
x  4> g  v;  O' will  be  used, 

7^e  [0.9517,.  0.9758]  ,  a  .ninth  order  autoregressive  model  with* 

0^  vt  will  be  used, 

7 ^ e  [0.9759,  Q v 9919] a  tenth,  order  autoregressive  model  with 

vi  0 '.will  be  used, 

7,^v  [0.9920,  1.00]:,  an.eleventfr  order  autoregressive  model  with 

v  0  will  be  used. 

Thus,  we  have  chosen  a  model  for  azimuth,  error  process,  y r),  is 
used  to  construct  elevation  error,  if 

7 2<r  [0,  0.6048],  a  third  order  autoregressive  model  will  be 
used, 

72e  [0.6049,  0.7580];,  a  fourth  order  autoregressive  model  with 

vi  0  will  be  used, 

72e  [0.7581,  0.8467],  a-  fifth  order  autoregressive  model  with 

<t> ^  •/  0  will  be  used, 

72‘.  [0.8468,  0.9192],  a  sixth  order  autoregressive  model  with 

/.  0  will  be  used', 

■6 

72e  [6.9193,  0.9676]1,  a  seventh  order  autoregressive  model 

with  je  0  will  be  used, 
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7„c  .[0/9677-,  0 .9757]',  an  .eighth,  order  autoregressive  model  with 

i>g  *  O1  will  be  used, 

[0.9758,  0.9838}\j  a  ninth  order  autoregressive  model  with 
Oy,  *  0  will  be  used, 

7^e  (0.9839,  0 i9919J ,  a  tenth  order  autoregressive  model  with, 

$10  ?-  0  will  be  used,  . 

7 2?  [0.9920,  1 .0] ,  an  eleventh  order  autoregressive  model  with 
<t> ^ ^  *  0  will  be  used. 


b)  Use;  normal  random  number  generator  to  generate  the  required 

2 

parameters'  p,  o's.,.  and  a  . 

c)  Using  a  polynomial  root  solver,  check  the  roots  of  Xp  -  o^Xp"''L 
....  -  $  v  -  0.  If  any  of  the  roots- is,  greater  than  or  equal  tU  1,  discard 
this  set  of  o's  and  select  another  group  of  parameters. 

d)  Let  X  denote  the  -azimuth  error  .process  and  denote  the 

elevation  error  process.  Then  according  to  the  models  and  parameters 
chosen  by  steps  a)  and  b),  we  can  simulate  X  and  Y  consecutively 

b  2^* 

by  generating,  normal  random  derivates-  a  from  N(0,  o  )  . 

t  cl 


•e)  If  the  perfect  aim  of  a  gunner  at  time  t  is  (A-  Ej.),  then 

our  simula  ted  coordinate  of  a  gunner  at  time  t  is  (At  i  V  ,  Efc  +  Y  ) . 
A  simulation  example:  ‘  *  * 
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Appendix  A. 

SIMULATION  PROGRAM  LISTING 
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Preceding  page  blank 


OGRAM  HA IN  7m/7<.  0PT*1 


FTN  6.3+P353  37/1G/73 


DIMENSION  ALIHIT  ( 111)  ,ELIHI  T  (10 1  ,XU  A(U>  ,XUE  (11 )  ,XLA(U)  ,  XLE(11 )  , 

- A-PH(tl>,PH£(ll*TSPm)  t3PE(ll)  >A  (98»W*&(  500) ,AZ15II)  ,CL(58i>- 

01  HENS ION  SPl(ll) »SBE1  111) .PHI  111) ,PHE(11) 

DIMENSION  IK(50i)) 

- -»HFA  •8/ll*wrt/ —  - - - -  - 

DATA  AUHIT/O. ,.6291 .. 6063, . 8549, . 9114 , .9636, .9517. .9759, . 992C ,1. / 

READ  1«1, (PM(I), 1=1,11) 

- -RE  A9  1S-I  rtSP«T-rI*i,lt-» -  - - 

REAO  U1, (PNE(I) ,1=1,11) 


- RrAO  tTt,+SPcttt»-I*t-,tl) - - - 

RE  AO  lul,  AHEAH  ,AHV,cHE  AN,£1  V,  VARA,  VVARA,VARE,  WARE 

- NPTS*~3flC- - - - - 

JUMP A  =  3 

C  OETERHINE  ORDER  OF  TIME  SERIES 

STtRT  *  0. 

...  - GALL  -RANN*M(i,>START-,RNA»-  -  ■- -  - 

CAU  RANNU.Nll.  ,STERT,RNE) 

DO  1 0  I  ■■  2,9  —  •  —  —  -•  ■  .... 

IF  I  (  RNA.GE  .ALIHIT  ( I)  )  .  AND.  (  RNA.LT*  ALIHIT  (  1*1) ).)  JUHPA=If  Z 

. IF((*Nc3&E-»EtIMIT<m.*N(hN-RNeTt-T.€UNmi'*-lTH--UtfMPe»lH* -  — 

13  CONTINUE 

fr  -  GENERATE  HU, PHI-S.WAR(A)  -  -  •  . . 

DO  332  I  =  1,11 

. SPKI1  ••  -50RHSPUH  - 

3m 2  S°CI(I)  *  SORT  (SPE (I ) ) 

..  ..  .  — 00--2v 

XUA(I)  3  p«ci)'  »  5,»SP1(I) 

XLA(I>-  »  Prl(I)  -  0,*  SPl  (IT  .... 

XUE<  I )  s  PME(I>  ♦  i.’SPEKt) 

20-XLETIT  =  PrtEfl)  -  5.  *SPE1H )  .... 

62  00  5m  I  =  1,11 

•  ■  ■  -»v-PHim •  «  g.o  -  -  .  -  -  — . .  . 

DO  m3  I  =  1,3 

— 65  CALL  NCWHAuJPHITIlMXLA  +  n.XUATD.PmiLfSPl-d))  - 
IF  { JJMPA.  ECl.C)  GO  TO  50 
•  ■  I  *•  JUHPA 


CALL  NORHAu<PHI<I1mXLA<I),XJ'TI)  ,PH{I),S“1(I)) 
- HN-«  -UOHPA - -  -  ••  - . .  •• 


GO  TO  »5 
-  53  MN  «  3  • 

55  R(l)  =  1.0 

90—57  I  -■  *1tNN 
57  R< !♦!)  3  -PHI ( I) 

. . -OAL-t— ROOTHHRHTtNN,  RR-tR  £ )  -  -  - - -  —  - 

IN0S4  =  i) 

00  5s  1  «  i.»NN 
IF  (ASS (RR ( I) ) • GE. 1 .)  IN0SA31 
56-  CONTINUE- 

IF (INOSA. He.,  1)  GO  TD  61 

.  ...  -pn*nT— 1j2-  -  • —  -  -  — . . . 

GO  TO  62 

-61  00-63  I  •*  1-,  11 
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OG'RAM  IAIN 


.fTN  6.0+P3S3 


07/1C/73 


-  — 63-PH€<B~»-6. - - - - — - -• 

00  65  I  =  1,3 

65  GAtt  NeRNM-WteXI-)  ,Xt*a+-rX0t«+-YPNe«-H5P€*fW4- 
IFUUHPE.tO.O)  GO  TO  66 

CAU.  NQRMAUPHETI)  ,XLETI),<JEtI)  ,PME ( I) ,SPEl < I) ) 

. _  _  ,  ii  k.  —  inu^r  .  .  -  -  ,  _ 

nn  —  9vn  tu 

GO  TO  67 

-  —66  -NH~»--3 -  -  - -  - 

67  R( 1)  s  1.0 

- - 90-62- --I  -  - - 

12  R(IH)  =  -PHc(I) 

- — — ROOT  S1R?8,N6|RR,RI)  : - - — . .  -•  —  — 

INOS;  =  a 

-  *  00- 66“ — t,WH‘  —  —  ...  . - ....  - 

IF<AaS(RR(I>>  .GE.1.»  ,IN0SE*1 

—6ft  CUNT-INUE-  ■  ~  . .  --  — - 

IFlINOSi.Ni.l)  GO  TO  69 

GO  TO  61 

69--AMV1  *  SORT (AMV)— -  -  •-  - 
EM  VI  =  SORT  TEMV) 

-  WARM  -*-SOR-HVVARAI  —  . 

VVAREl  =  SORT  ( VVAPEJ 

- Xfc-HA— »  -AM€AR-^r»*HVlr  -  - - - -  - 

XLMF.  =  £Mc«N-3.»tMVl 

-  XU*A  -*~AM£AN+3.*AMVi  -  ■  —  —  • 

XUME  •  EMb\N4  3.»EHVl 

GALL  NLAHftL'TAMN,XLHA,XyHA,  ANtAN,  AM¥t>-  -  - 

CALL  NCRMAL(EMN,XLM£  , XUME ,  E MEAN',  EH  VI) 

XLVA  ■*— VARA-  —  -5«*WARA1  -  - . . . 

XL VE  =  VAR.-.  -  5. *VVARt  1 

XU VA  a  VARA  ♦  5.»VVARAl 

XUVE  =  VARt  ♦  5.*VVARE1- 

GALL  huRRALR AVAR  , XLVA,  XDVA, VARA,  VV ARAD 

CALL  NORMAl(EVAR,XLV£,XUVE,  VARE, VVAREl) 

-C  ~G€HERAlt--MT+M*.  -  . .  —  -  - 

AVAR1  =  SORT ( AVAR) 

EVAR1  =  SORT  (EATAR)  -  ~ 

XLAA  =  «5 « »AVAR1' 

XL  At  =  -?,*EVAR1»  •  — 

XU  AA  =  5, *AVAR1 

-  — .  -XUAs  --5.»«VAR1--  - 

00  8J  I  =  l.NPTS] 

GAL-L  NGRMtfLX AR  I)i,XL4A,XUA4-, } .,  AVARD- —  •  - 
83  CALL  N0RHAL(AEU!iiXLAE,XU4E,S.»EVARD 
C  PRIM  I.W-IALU  SELECTED  VALUES  •  - 

PRINT  1 

•  -  .  -  4F-WUHRAx^4r41-.,0  40  GO-  .  .  . . - 

PRINT  2  ,  AMN,AVAR,  ( PH  1(1). ,1  =  1 ,31  ,  JUMPA.PHI  (JUMP A) 

GO  TO-  13S .  ■  -  — . 

93  PRINT  3, AMN, AV/AR,  (PNI(I)  ,1  =  1,3) 

—  -TrOv- -IF.WUKRE-,*0«l  CO  TO.  US--  -  - - 

PRINT  t.EMN.c.TAR,  (PrtETI)  *1-1,3)  ,  JUMPE.PHI  ( JUMPS) 

- G0-TQ-I2S-  - - - - ~ 

itu  PRINT  5,ErtN,EVAR,TPHE(U,IU,3) 

-  - 5 -WRMAT.(3XJi/.WVATI0N*F6.iH‘4:r5.A)  -  -  -  -  . 
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0J3»H  MAIN 


FTN  6.0»°353 


i37/;:0/7J 


-tH i  PSINT-o - - 

COMPUTE  AZ  ANO  EL  VALUES 

- AZ  <!)■  =r-*MN  ♦  AI1> - -  - - — - 

AZ<2)  =  ANNAPHI  (1)  *(  AZ  (D'ANN)  »A(2) 

—  ArZ~tOt~~~-A>tWtPM-I,<l't'+tA'Z-<Eh^AHNt'*^Hi't?t~*'tAt'<  1>»AHMT  tMJt - ~ 

IF ( JUHPA.  E 0.0 )  GO  TO  300 

00  -3 rt-'-t-tr , -JtWP A - -  -  - 

3  if  AZ  (I)  *  AHNfPHIUI»(AZ<I-il-AMNI*PHI(2»MAZ(I-2‘-AM'0+PHH3** 

. A  fAZ(I-3T-AW0— t-AH>-  — - - ; - — 

J  J.  JUMP  A  ♦  1 

00-  330--I  *  J.NPTS  - - - - 

3  30  AZ(I)  a  AMNt-PHU  1)  *<  AZ  ( I-l)  -ANNI  *PHI  (2T*l  AZ  tI-2l“AMV)  tPMI  (3)  * 

-  »-  < *g**-3h-*MHt+PH-fi-3UMP»»-H-AZ-<  I-UUHPAi  -ANN >■■•■♦■  A  ,W - - 

GO  TO  3<*J 

3jC  00  320  I-  »-<*,NPfS - - — - - -  - 

320  AZ  ( I )  =  AMN+PHI (1T*(AZ(I*1)  •  AHNJ ♦PHI  (2) *( AZ (1-2) -AMM)  ♦P-II  C 3) * 
- A-<AiH-i-3»~AtW»  HH-tt -  --  - - - — 


-jii-ftm  -—ehn+aeh* — - -  -  —  - - - 

EL  (2)  =  EHN»PHEU)*(,ELfl)-E'1N)  *A£  (2) 

=L«»;  »  6NN*PM£UP*I-£t-I2)-EmW*PH€-(2T-M€L-(-l)-EI«H-»AE(3>-  - 
IF  ( JUMPS.  EQ.U  GO  TO  SCO 

00  <tlv  I  =  S.  JUMPS- . -  ■  -  - -  - - 

,10  EL  ( I)  =  £MN»PrtE(i}*<eL(I-il  -£HN)  *PHE (2)  *  (EL  (I -2) -END  ♦PMS  C  31  * 

A  '6eLU-W”cHN>**€fW  —  “  -  - - - 

J  a  JU.4PE  ♦  l 

00  -I  *  J.NPTS  . -  -  --  - 

♦  hC  EL(II  =  £hn«-PH£(1}*(EL(I”I)  -EMN)  «-PH£  (2>*(EL(I  — 2>  —  £KXi  ♦•?■(£  C  3)  • 

A  <EL<I-3)-EMN)+PHEU;>NP€>M£L-<I-JUHP£)-EHN)*AtfI)~  -  ■ 

GO  TO  -.33 

St,}  OO-S-26-  I— SrNPT-S-  - -  - - —  -  — -  -  -  - 

420  £L  ( I  >  a  £MN*PH£  (1)  *4  EL  (I-l)  -EMN)  .PHE  (Z)MELCI-Z)-ENM)  ♦PrfE  I  3)  * 

A  <ELU-3)-£HN)  fAEd)-  —  .... 

PR  I ‘41  VALUES 

45C  00  S31  I  *  l.NPTS  ... 

433  IK  ( I )  =  I 

-  13  »-  0  - -  -  . - 

IE  -=  C 

IFINPTS.LT.229>  -GO-TO  S51  --  .... 

13  =  Ii  ♦  i 
It  *  Id  +  SS- 
•» 5*.  00  *.52  I  »  18.  IE 

-.52  PRINT  (M*« t-rA24-I4rCL-MK-I»HI-»35 >-rA2-(H4S ) rEL-H ■*•»»> rK4 1 ♦!*»»*- 
i  AZ  II  HH>  tctUMUl!  IK  11*165  )»AZ  (1*165)  ,EL< 1*165) 

PRINT  i SS  * 

PRINT  o 

IF  <  <NP  F5-«»*Ic)  *LT  .2201  GO  TO  S51 
10  =  <*•!£  ♦  i 

££-.  HW— *  ?<i  - -  - -  -  -  — - 

GO  TO  *.5*. 

.01  NlEP-T«nPTS-S*I£ 

NCOLsNlEFT/Ss 
NEK  •  NLE F  T~NG0L-*»5 
13  =  S  MS  ♦  l 

IE  *-  id  ♦  -is  —  •  —  * - -■ 

GO  TO  (<.55.«,S6. .57), NCOL 
.55  IF  (NEX  »£Q.  o>  SO  TO  -.61  -  - 
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353AM  MAIM  7h/T+  0PT  =  1  FTH  A,a+?353  07/l«/73 


+IEX  —  1-3 -+ -Ne* - 1-  -  -  -•  - - - 

00  *58  I  =  I3,NcX 

♦  6«  PRINT  (bIKUHA2<I»,ELMI»,KU+55T  ,AZTI*55),et<-I+5$T  - - 

NXi  =  NcX  »  1 

06-<r59'I— •  M*liK .  --  -  -  •  —  - 

*03  PRINT  3,IK(IJ  ,AZ(I>,ELtI) 

*  5O-"T,9-l-0fl'8-t,  ■ — - -  -  —  — —  — ■ —  ■  *  —  —  — , - - — 

-,61  00  ao2  I  ;  18,  It 

♦62-  PRINT  »,I*«U  ,AZ«>reU-a>  —  —  •  •  -  -  ------ 

GO  TO  lu/,. 

♦-53  If  THtXftUml  50,  -+31  ■  -  ■-  —  -  -  — 

N-X  =  10  ♦  NiX  ■  1 

•00 — tov- I.  -* — Ttfr-R  -  --  - -  -  -  *  -  --  -•  —  -  —  — 

■,6-,  PRINT  3,  IX  (I)  ,JkZ<I>,ELfI),l<(I»55»  ,AZTU55>  ,£  Li  1  +  55)  ,!<(!  + 113)  , 

1  AZa>il«»,efc-«+ll3) 

NXI  =  NiX  ♦  i 
00-  to5  I  =  4X1, IE 

-.65  PRINT  8,*KU>,AZUI,ELII>,I<<I*i5)  ,AZ  ( I+55>  , ELI  1+55) 

GO— TO  --I d -)-♦ v  — -  —  -  --  -  -  — 

♦  63  00  <,56  I  s  18,  It 

♦  65  PRINT  9,IX(If,  AZ<H,€LM>,I<lI*»5t  ,AZII  +  »5I  ,tUI*55» 

go.  to 

♦  67  IMNeX.EO.j)  GO  TO  ^73 

N;X  =  i<3  ,  NtX  -  1 

00  ♦-?♦  I  «  18, HEX  ■—  - 

-,?♦  print  3,  IK  1 1)  ,  AZ(I»,tL<I),I<(I*55)  ,AZ  <1+55*  ,£L(  1+55) ,  I<  <  I  ♦  110 )  ♦ 

1  AZi  1  +  111  >  rt-LT  I+Uv,)-",lK<-m65>  ,AZTI+-165»  »€li  1  +  165) 

NX  1  =  N ;<  »  l 
00  A/-,  l  =  NXI, IE 

,76  PRINT  0,iKlI),AZU),tLU>,I<U+55> ,AZ  U  +  55)  i£L  U+55)  ♦  KU  +  IU  ) , 

1  AZTI  +  H'i), £1(1  +  1 13) 

GO  10  *j.  ,< 

,73  00  ,76  I  =  13,  It 

,76  POINT  d,I«-!),AZ(ll,tL<I>,t<<I*55>  ,AZ<  1+35)  ,EL  U+55) ,  I<  ( I  ♦  110  >  , 
l  AZ(I*ilvl«cL(I  +  113> 

1  FORMAT ( lHl. 2s X* INITIALLY  SiLcCTiO  VAlJES*/13X*MEAN*2X*AVAR*2X 
4»PHI-*»tX»PNl2»2X»PHT3»2X»JONP»2X*PHI-TOWM0)  »> - 

2  FORMAT  (3X»AZI4UTH  *F6. A, 4F6, A , 3X12, aXFo. A) 

3  FORMAT  (3X  +  AZIMUTH  *F6,At^FS,A) 

♦  FORMAT (3X*tLtVATI0N*F6.4,,P5.A,3xi2,AXF6,+) 

3  FORMAT (//,3X* TINE ♦3X*AZINJT M*3X»-LEVAT ION*, 31 AX* TIN;* 3X+ AZINUTH+3X 
1  +SL£V«TI0n*l) 

7  FORMA T-<3XI4r6-XF6rA-T*XF6. AT . .  -  -  - 

9  FORMAT (♦XIG,1XF9.A,5XF8.a, 3 (5X13, IXF3,A,5XF8.A>  > 
lvl  FOP.NAT  (5XHF6.U 

iU  FORMAT  (/// ,  1. X* INOSA  =  1,  RtITSRATEO*) 

,u3  FORMAT  l///»l(,X*INOSi  =  l ,  R; I TERATtO* > 

It,  FORMAT (1M1 ) 

I'OuOi-  ST-OP  -  -• 

ENO 
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IN:  NONNA L  74/7,  0PT*1 


FTM  4.9*P3$3  07/1C/73 


■  ~SO«WtfriNe-Ni«N*l-(  RX»XLrXTMHt»ttGNAT - 

SUBROUTINE  NORMAL (  *X,  XL.  XU.  NU,  SIGMA  » 

■  'OC9CRh*HON - 


THIS  ROUTINE  PRODUCES  NORMALIZED  RANDOM  NUMBERS  NITNIN 
-A-iHN»i-^R«IPK»-  Br-TNE-UBiRt - 


INPUT 

-  • — -i — XL- - -tONER^rlNIf— 0P-N9RNAL  -ONRVE - - - 

2  XU  UPPER  LIMIT  OF  NORMAL  CURVE 

«.  SIC  STAN3ARO  DEVIATION  OF  NORMAL  CURVE 

OUTPUT - 

1  RX  THE  NORMALIZED  RANDOM  NUMBER 

01 HENS  ION  Fill 261 

■  REAL— HU—  - - - - 


DATA  0 1.0  T1.NTP.NPT.SQ2PI/. 01. .34. 126. C  «.  396942263/ 

- SL?,"»~Nd - 5r9~*  SIGMA-  -  - - 

SUE  s  HU  t  5.0  *  SIINA 

IfH  -XL  -ri.fn~ SU5-  vANOr— XU  .  3Ei-  9t5  I — SO-f-0  H —  ~  ■ 

PRINT  *J1,  SLS •  SU5,  XL.  XJ.  MU.  SIGMA 

1  »  ANU’FIC.2  •  —  PROGRAM  TERMINATED.*///*  XL  «*F1J.2.5X 

—2  —  *-*U-**-F10.2rW*-Htf- »*F19-r2>3X*-SI6NA-«»F*6» 

STOP  lc3 

li-  IF<NPT.-Ne.i»-SU  TO -2  - --  - -  - 

T2=0T1 

FZ  III  *  .5"  •  - -  - . .  . . 

FC-FZ 1 1) 

Tan.  - —  .  • — * —  - ■  — 

FP=S02HI*EXP(-.5*T*TI 

•OT2*9T*.5-  *  ■-  '  —  •*  -  —  — 

NPW1 

1  CONTINUE  -  -  —  - 

T  *  WOT 

F  «  SII2PI *EXPI •»5*T*TT-  - 

FC=FC  ♦  0T2*IF»FPI 

FP=F  •-  •  -  •  -  ■•  ... 

IF(A3SlT-T2I.GT..C5ilt  GO  TO  1 
NPWNP1  t  1 
FZINPTM'FC 
T2M 2*0T1 

IF(NPT.LT.NTP)  GO  TO  1 

2  CONTINUE 

SH  »  RANFIRI  ■  ■  - 

W9.N 

IFIRN,lT..»)  R*l.-RN  . 

IF  (R.G  TiFZ  INPT  1 1  GO  TO  2 

00-3-Hl.NPT  --  —  ■  -  - - 

IFIR.GT.FZUH  GO  TO  3 
IFU.EQ.lt  30  TO  **  •• 

X*I-I 
X*X*OT 1 
Xl»X-OTl 

RX  *  Xi-*-(R-FMI'lll*BTt/T«2iI»-fZU-ltt- 
GO  TO  i 

3  CONTINUE  ■  •  • 


II.NE  NORMAL 


7,/7  •*  0°T =1 


FTN  ..TtPSSS 


PRINT.-lJ* 

IvQ  FORMAT  U3Ho  ERRgR  0ISRN1 
STOP 

4  30NTINUE 
•RX»E. 

5  CONTINUE 

■  -  IFIRN.LWtW-RX«-RX- 
«X»  SIGMA*RX*MU 
•IFIRX.GT.XUt  GO  TO  l 
IF(RX.LT.XL)  GO  TO  2 
RETURN  ■ 

ENU 


~ 

NORM 

1C 

--NORM 

•26— 

NORM 

30 

■  NORM 

-49 

NORM 

5C 

NORM 

70 

'ilftny 

WIMI 

se- 

NORM 

90 

AA  (>f>M  4  n  A  _ 

NORM 

110 

-NORM 

126- 

NORM 

13D 

NORM 

146 

NORM 

ISO 

NORM 

166 

NORM 

17C 

NORM  IMG  - 

NORM 

lEC; 

MORtf  2V& 

M0RM‘Sjl(! 

NORM 

i?X 

NORM 

23C 

NORN  249- 

NORN 

25C 

NORN 

26C 

NORN  271 

NORN 

289 

NORN 

290 

NORN 

300 

NORM 

31t 

NORN 

323 

NORM 

33C 

NORM- 

340 

NORM 

300 

-NORM- 

369 

NORM 

370 

NORM  385 

NORM 

39C 

NORM 

400 

NORM 

410 

NORN-426 

NORM 

430 

NORM 

449 

07/K/73 

-NORM 

450 

NORM 

460 

NORM 

479 

NORM 

48C 

-NORM- 

490- 

NOR** 

S3C 

-NORM 

519 

NORH 

520 

-NORM 

•530 

NORH 

54C 

NORM 

•559 

tiORM  560 
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